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We investigate the possible effect of the residual annihilations of heavy particles after freeze-out 
from equilibrium in the early universe. An error in previous work on this subject is pointed out 
and the correct method of solving the Boltzmann equation for this case is developed. For Majorana 
particles there are significant differences relative to previous work. 
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&\ , Introduction: In the early universe particles are subject to both pair-annihilation and elastic scattering interactions. 
At sufficiently early times these interactions are usually strong enough to keep the particles in complete thermodynamic 
equilibrium. However, as the universe expands and cools, interactions become weaker and eventually any particle 
species will decouple from equilibrium. A quick estimate of the freeze-out epoch can be found by comparing the 
expansion timescale of the universe, -ff -1 , with the interaction timescale for the given species, [Q. If the interaction 
■ timescale become longer than the expansion timescale the species decouples so that the freeze-out criterion is 

c> (i) 



For massless particles with standard weak interactions one finds for instance a freeze-out temperature of Tp ~ 1 



> 

MeV. If the species is absolutely stable against decay, its present contribution to the energy density of the universe 
£ — ■ can basically be found from the ratio m/Tp (except for the possible dilution of the species by entropy production 
\ subsequent to freeze-out). This can be used to constrain models with stable massive particles by demanding that 
their present day energy density does not overdose the universe The standard example of this is the Weakly 
Interacting Massive Particle (WIMP), a hypothetical heavy particle species (for instance the lightest supersymmetric 
particle) which could make up the dark matter of our universe. 
Oh| Even if a given particle species violates this constraint it may still be allowed if it decays on a sufficiently short 
timescale. There are, however, other possible ways of constraining such scenarios. For example, a massive decaying 
neutrino will in general change the outcome of Big Bang nucleosynthesis (BBN) ||. Further, if its decay products 
c/2 ■ have electromagnetic content the model can be ruled out by other observations [H]. 

Now, even for a stable particle there may be other possible constraints than the energy density argument. After 
freeze-out there continues to be some residual annihilations of these particles. If their annihilation products interact 
k>( \ electromagnctically and are sufficiently energetic they may photodisintegrate light elements and ruin the agreement 
between BBN theory and observations. For this reason BBN can be used to constrain the interaction strength and 
mass of WIMPs. This has been done several times in the literature, both for high mass particles (to > lGeV) [0-0], and 
for more moderate masses [0j8|. In all cases, the so-called integrated Boltzmann equation has been used to determine 
the annihilation rate of heavy particles. This equation assumes that all particles are in kinetic equilibrium, but, as will 
be pointed out, it has not been taken properly into account that any massive particle species which decouples from 
chemical equilibrium eventually also decouples from kinetic equilibrium. Using the standard integrated Boltzmann 
equation thus leads to wrong estimates of the annihilation rate. 

The integrated Boltzmann equation: The approach to take when calculating the freeze-out of massive species is to 
use the Boltzmann equation 



df = C a „n + C eU (2) 

where 

»t 9 f IT d f 

d! = Tt- Hv W (3) 

H being the Hubble parameter H = R/R. On the right-hand side, C ann and C e \ represent annihilations and elastic 
scatterings respectively. If one assumes that the interaction is CP conserving these collision terms can be written 
generically as [[!) 
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C coU [/] = —J d 3 p 2 d 3 p 3 d 3 p4Hflj2,f3j4) 

xSJ2 Mh^ 4 (Pi+P2 -ps- Pi)(^)\ (4) 

as long as we are only concerned with 2-body collisions of the type 1 + 2^3 + 4. Here we have d 3 p — d 3 p/ ((2ir) 3 2E). 
S is a symmetrisation factor of 1/2! for each pair of identical particles in initial or final states jnj, and ^2 \M\ 2 is the 
interaction matrix element squared and spin-summed, pi is the four-momentum of particle i and A(/i, f 2 , f 3 , f&) = 
73/4(1 ± /i)(l ± f 2 ) — j "1/2(1 ± /s)(l ± /i) is the phase-space factor with + corresponding to Fermi-Dirac statistics 
and — to Bose-Einstein statistics. 

If one then assumes scattering equilibrium for all particles involved the elastic scattering terms all equate to 
zero. Further, a significant simplification is possible if one uses Boltzmann statistics for all particles instead of the 
appropriate quantum statistics, so that 

1 ± / - 1, (5) 

where, again, + corresponds to Fermi-Dirac statistics and — to Bose-Einstein statistics. This reduces the phase-space 
factor to 



A(/l,/ 2 ,/3,/4) = / 3 /4-/l/2 (6) 

The last assumption is that the annihilation products are in full thermodynamic equilibrium so that 

hh = e-^+ E ^' T = e- E ^ T e- E ^ T = ./,,„/,.,„ (7) 

By integrating over momentum space for the incoming particle one then arrives at the well known integrated Boltzmann 
equation |l|U 

n = -3Hn- (av)(n 2 -n^ q ), (8) 



where (av) is the velocity-averaged annihilation cross-section, defined as |l[|llf 

(av) = -\- I d 3 p 1 d 3 p 2 d 3 p 3 d 3 p 4 f 1 , cq f 2 :Cq 

xSY, M 2 i2^ 4 (Pi+P2 - P3 - P4)(2tt) 4 , (9) 
and n eq is the equilibrium number density 

" c q = 0~^2 / P 2 fcq{p)dp. (10) 

Below we list the assumptions going into the integrated Boltzmann equation, Eq. (g): 

1: The heavy particle species is in scattering equilibrium 

2: All annihilation products are in thermodynamic equilibrium 

3: Boltzmann statistics instead of Fermi-Dirac or Bose-Einstein statistics 

4: CP invariance. 



At very high temperatures the heavy species will be in complete thermal equilibrium. However, as T drops below 
the rest mass m, the rate for pair annihilation and production becomes exponentially suppressed, V cx e~ m l T , quickly 
leading to freeze-out of chemical equilibrium at some temperature, Tp A . 

However, the rate for scattering on other relativistic species is not exponentially suppressed. This means that kinetic 
equilibration is maintained long after freeze-out of chemical equilibrium, but at some point there will be a freeze-out 
of scattering equilibrium too. From that point on the distribution will remain fixed in comoving momentum space 

f(p,T) = f(p F T Fs /T > T Fs ), (11) 

where Tp s indicates the temperature for freeze-out of scattering equilibrium. At the point where scattering freezes 
out the distribution has the form 
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f{p,T Fs ) = e -(m+P /2m-n)/T Fs _ ^ 

At later times, this distribution is no longer of equilibrium shape. Rather it is an equilibrium distribution corresponding 
to a new temperature parameter 

T.=2>£. (13) 

1 F 

Since the integrated Boltzmann equation only applies to species in kinetic equilibrium it is not apparent that it should 
be applicable to a treatment of post freeze-out annihilation. Note, however, that for a non-relativistic species after 
freeze-out, / cq <C /, so that pair-production can be neglected. Thus, the annihilation term in the full Boltzmann 
equation reads 



xSY, Mi2^34<5 4 bi +P2-P3- Pi){2K)\ (14) 

where / is a distribution function corresponding to scattering equilibrium with a temperature T*. Thus, one can still 
use the integrated Boltzmann equation in the shape 

h = —3Hn — (av)n 2 , (15) 

but with a different definition of (av). An annihilation cross section which in kinetic equilibrium is proportional 
to T n becomes instead (av) oc T™. In general we can therefore state that after chemical freeze-out there will be a 
temperature region, Tp A >T> Tp s , where the equilibrium assumption (av) oc T n applies. Then at T = Tp s kinetic 
equilibrium decouples and the thermally averaged cross section follows the new relation (av) oc T™. 

All previous work has used the equilibrium assumption, (av) oc T", even after kinetic freeze-out, which gives a 
significantly different temperature dependence of the annihilation rate. For Dirac particles, n = 0, and the result is the 
equivalent to the standard one. For Majorana particles, n = 1 [], and there is now a big difference between (av) oc T 
and (av) oc T*. The cross section falls off much faster with decreasing temperature so that residual annihilations are 
weaker than found in the standard calculation. 

How big this effect is depends on the difference between the kinetic and chemical equilibration rates. In general 
scattering and pair processes arise from the same coupling, so that the fundamental interaction strength is the same 
for the two different processes. The thermally averaged scattering cross section will therefore be (av) oc toT. However, 
what we are really after is the thermalisation rate for the heavy particle, not its interaction rate. In general a non- 
relativistic particle scattering on a massless species only gains or looses a fraction T/m of its energy. Therefore the 
real thermalisation cross section should instead be (av) oc T 2 . Since the relativistic species have number densities 
going as T 3 the overall thermalisation rate can be estimated as 

T T oc T 5 . (16) 

This rate is in general much larger than the annihilation rate because it does not have an exponential suppression 
factor in the number density of scatterers. The ratio of annihilation to thermalisation can be written approximately 
as 

^H?f V ^, (17) 

as long as scattering equilibrium still holds. For a Dirac particle decoupling from chemical equilibrium when m/T ~ 20 
this ratio is roughly 10 -4 at chemical freeze-out. 

Applications: There are a number of cases where the above formalism may be applicable. The standard example, 
of course, is the massive neutrino. Frieman et al. H have developed an analytic estimate of the effect of residual 
annihilations Their results only apply to masses below the threshold for neutron and proton production (to < 1 



1 This applies as long as the massive species is non-relativistic and the annihilation products are massless. 
2 Their treatment of the photon cascade processes is not numerically accurate and lacks several relevant processes, as pointed 
out in Ref. Q . However, since we are not after high precision in our calculations we shall use the results of Ref. M anyway. 
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GeV), but will serve well to illustrate our results (for a treatment of higher masses, see for instance Refs. HH). 
Energetic particles injected into the cosmic plasma after BBN can photodisintegrate the light elements. The primary 
effect is to fission 4 He to produce 3 He and D. Since the primordial element abundances are quite well determined one 
cannot, for instance, tinker too much with the ratio of helium isotopes without coming into conflict with observations. 
The primordial value of 4 He has been determined by Olive and Steigman |D3] to be 

Y P = 0.232 ± 0.003(stat) ± 0.005(syst), (18) 

where Yp is the mass fraction of helium. Other determinations have been slightly different, but since we are not 
interested in great accuracy we shall just use the above value without further discussion. For 3 He the primordial value 
is bounded from below by observations of the local interstellar medium to be |l4| ] 

N(H) ~ N(H) ~ V ' 



One then has an upper limit to the ratio 3 He/ 4 He of 



* 2 x 1(r3 ^ 

As an example we take the annihilation of very massive neutrinos with standard weak interactions. For both Dirac 
and Majorana neutrinos one finds that for a mass of the order 1 GeV, the chemical freeze-out happens at xp ~ 17 giving 
a chemical decoupling temperature, Tf a , of 60 MeV. For Majorana neutrinos the scattering equilibrium decoupling 
happens at 7> s ~ 1 MeV. 

Using the constraint above on the helium isotope ratios Frieman et al. |?J have derived an estimate of how the mass 
and interaction rate of a given particle can be constrained 



TOMcV 



> 1.0 x lQ- 2 ex 2 , 



2-1 !/("+!) 



CO 

2.6 x 10 4 5 ^ 1/(n+1) 



g*(T F )n B h 2 (n + l)(n+D 



(21) 



Here, the following parametrisations have been used: (av) = ao(T/m) n and the number of massless degrees of freedom 
is defined as = ^jT^ 4 pp. B is the branching ratio to electromagnetic annihilation products and Xf = m/Tp. 

The above constraint was derived assuming that the equilibrium version of the Boltzmann equation is correct. 
However, it turns out that we can still use it for Majorana neutrinos. The photo destruction of Helium only begins 
after T = 1 MeV § which means that at the time of 4 He destruction the Majorana neutrinos have also decoupled 
from scattering equilibrium. This means that (av) oc T 2 . However, we cannot just use the above relation taking n = 2 
because that would assume that the scattering equilibrium decoupled at the same time as the chemical equilibrium. If 
we first rescale the overall annihilation cross section with a factor -0 = Tp A /Tp s ~ 60 then we can use the constraint, 
Eq. (^l|), by using the substitutions 

cto -> coXFip, (22) 
n = 1 -> n = 2. (23) 

Plugging in numbers one finds that for a Majorana neutrino with standard weak interactions there is a lower limit 
to the mass of 

m v > 150MeV (n = 2), (24) 
which should be compared to the value one finds by using n = 1 

m v > 300MeV (n = 1). (25) 
For a Dirac neutrino the mass limit is of course unchanged and is found to be [Q 

m v > lGeV (n = 0). (26) 
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The above limits have been derived using fi^/i 2 = 0.05, B = 0.5. 

Thus, using the correct version of the Boltzmann equation weakens the mass limit on heavy Majorana neutrinos by 
about a factor of two. In itself, of course, the above limits are not that interesting because massive neutrinos heavier 
than m ~ 20 MeV are already ruled out if they couple to the weak interaction with normal strength. However, this 
weakening of the mass limit can be expected to apply generally because it only depends on the different temperature 
dependence of the annihilation cross section. 

Other groups have studied the effect of WIMP annihilations by means of numerical work |5|(| , taking into account 
for example the production of neutrons and protons before or during nucleosynthesis. This approach, however, makes 
it very difficult to estimate the effect of using the correct Boltzmann equation, though in general the impact of WIMP 
annihilations will, as mentioned, be less drastic than usually thought because the cross-section decreases much faster 
with decreasing temperature. This will make any bound on particle masses and/or interaction strengths looser than 
has previously been found. 
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